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Abstract 

Background: We are exploring the viability of a novel approach to cardiocyte contractility assessment based on 
biomechanical properties of the cardiac cells, energy conservation principles, and information content measures. 
We define our measure of cell contraction as being the distance between the shapes of the contracting cell, 
assessed by the minimum total energy of the domain deformation (warping) of one cell shape into another. To 
guarantee a meaningful vis-a-vis correspondence between the two shapes, we employ both a data fidelity term 
and a regularization term. The data fidelity term is based on nonlinear features of the shapes while the 
regularization term enforces the compatibility between the shape deformations and that of a hyper-elastic material. 

Results: We tested the proposed approach by assessing the contractile responses in isolated adult rat cardiocytes 
and contrasted these measurements against two different methods for contractility assessment in the literature. 
Our results show good qualitative and quantitative agreements with these methods as far as frequency, pacing, 
and overall behavior of the contractions are concerned. 

Conclusions: We hypothesize that the proposed methodology, once appropriately developed and customized, can 
provide a framework for computational cardiac cell biomechanics that can be used to integrate both theory and 
experiment. For example, besides giving a good assessment of contractile response of the cardiocyte, since the 
excitation process of the cell is a closed system, this methodology can be employed in an attempt to infer 
statistically significant model parameters for the constitutive equations of the cardiocytes. 



Background 

Introduction 

Cardiovascular research based on enzymatically disso- 
ciated cardiocytes has been fundamental for the discovery 
of the mechanisms that govern the heart. The use of the 
cardiocyte as the basis for cardiac functionality has pro- 
vided some of the most revealing information regarding 
heart function. Among the many findings, it has revealed 
the crucial molecular changes that occur during patholo- 
gical conditions of the heart. The details regarding the 
excitation-contraction coupling, calcium transient signal 
(movement of the calcium ion Ca 2+ ), gene and protein 
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expression, and contractility are all important mechan- 
isms and functions that can be readily studied in the iso- 
lated cardiocytes at all stages of development and they 
are routinely performed during research studies [1-6]. 

Contractility in adult cardiocytes is commonly inter- 
preted as the ability of the cardiac cell to generate force 
and to shorten. Some of the different methodologies 
devised to study the contractile process include laser dif- 
fraction [7], photodiode arrays [8], scanning ion conduc- 
tance microscopy [6], and those employing microscopic 
cell image analysis [9-12]. Historically the most widely 
used methods have been those involving cell image analy- 
sis, although all the methods show some positive and 
negative characteristics that are worthy of attention. 

Methods such as the scanning ion conductance micro- 
scopy require elaborate and expensive equipment [6]. 
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This technique, combined with laser confocal micro- 
scopy, is one of the few methods that has been capable 
of providing a measure of cardiocyte height during 
contraction. Other methods, such as light diffraction 
techniques, have been applied to the study of muscle 
mechanics since the nineteenth century with relatively 
high reliability. Nonetheless, they are very dependent 
upon several factors including the temporal resolution 
of the detection system and optical artifacts [2]. The sar- 
comere striation pattern analysis method has also been 
used as a way to quantify contractility. This technique 
can achieve high temporal resolution with the aid of 
charge-coupled device line array detectors and it pro- 
vides a measure of individual sarcomere lengths along 
the cell [13,14]. A drawback of this method is its vulner- 
ability to errors introduced by slight rotational and 
translational changes that normally occur during cell 
contraction [2]. 

Although the image analysis methods have been 
widely used with relative high reliability, the results they 
provide are often prone to the introduction of error due 
to the aforementioned rotation or vertical and horizon- 
tal displacement of the cardiocyte during contraction. 
Our proposed approach aims to provide a cardiocyte 
contraction analysis method that successfully captures 
the full extent of the contractile behavior, while mini- 
mizing the need for elaborate equipment and the effects 
that the cardiocyte's movements have on the acquired 
signal. 

Previous Work 

A widely used video-based method to measure contrac- 
tion in adult cardiocytes involves a device capable of 
capturing the extent and rate of length shortening 
between the cell's ends, the so-called edge detection 
method [15-17,11]. In this technique, a single video ras- 
ter line oriented along the longitudinal axis of the cell 
shows high contrast at the cell boundaries. These serve 
as tracking points and their separation distance corre- 
sponds to cell length. The method generally produces 
satisfactory results and has been a broadly used 
approach for measuring contractile responses of adult 
cardiocytes for over twenty years [2,16]. Some practical 
difficulties have been identified during the implementa- 
tion of the edge detection method for measuring adult 
cardiocyte contractility [2]. The changes in cardiocyte 
geometry, dynamic torquing, and rotation can lead to 
errors in the measurement [2,15,16,5]. Figure 1 shows 
frames extracted from videos of adult cardiocytes 
depicting contractions. 

In a previous communication, Bazan, Torres-Barba, 
Paolini and Blomgren [18] described a computational 
pipeline for the comprehensive assessment of contractile 
responses of enzymatically dissociated adult cardiac 




Figure 1 Frames from contracting cardiocytes. Frames extracted 
from videos of adult cardiocytes depicting contractions. These are 
typical rod-shaped cardiocytes isolated from an adult mammalian 
heart. 

\ J 



myocytes. The methodology comprises the following 
stages: digital video recording of the contracting cell, 
edge preserving total variation-based image smoothing, 
segmentation of the smoothed images, contour extrac- 
tion from the segmented images, shape representation 
by Fourier descriptors, and contractility assessment. The 
physiologic application of the methodology was evalu- 
ated by assessing the overall contraction in isolated 
adult rat cardiocytes. The results demonstrated the 
effectiveness of the approach in characterizing the more 
appropriate, two-dimensional, shortening in the contrac- 
tion process of adult cardiocytes. The authors in [18] 
compared the performance of their method to that of 
the aforementioned edge detection system. The method 
not only provided a more comprehensive assessment of 
the myocyte contraction process, but can potentially 
eliminate the historical concerns and sources of errors 
caused by myocyte rotation, bending, or translation dur- 
ing contraction [2,9,19]. 

In this paper, we are exploring the viability of a novel 
approach to cardiocyte contractility based on biomecha- 
nical properties of the cardiac cells, energy conservation 
principles, and information content measures. The pro- 
posed methodology was inspired by the works of Byung- 
Woo Hong et al [20-22], School of Computer Science 
and Engineering, Chung- Ang University, Seoul, Korea; 
and Andrew D. McCulloch et al [23-26], Department of 
Bioengineering, University of California San Diego, La 
Jolla, California. 

Methods 

Cardiocyte Shape Representation With Integral Kernels 

We wish to retrieve the information embedded in the 
shapes of a contracting cardiocyte like the ones shown 
in Figure 1. In other words, we want to analyze closed 
planar regions D c R 2 , and their boundaries (finite peri- 
meters), as depicted in Figure 2(a). We will describe 
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Figure 2 Cardiocyte shape representation, (a) Representation of 
the cardiocyte shape by a closed planar region, D. (b) Binary 
representation 5 of the cardiocyte shape. 





Figure 3 Cardiocyte shape features (asymmetric). Examples of 
the shape features R a {x\S) = S (x) {K a * (1 - 5 (x))) for (a) o = 15 
and (b) a = 25. 



these regions by binary images, Figure 2(b), composed 
with a suitable class of (continuous and invertible) 
image domain transformation. The binary images will be 
represented by a characteristic function 



S(x) = S D (x) 



if x e D 
if x & D, 



(1) 



defined for x e CI c R 2 , with D <= n, where CI is the 
rectangular image domain. 

We will define the multi-scale nonlinear features R G1 
as the convolution of the shape S with a family of ker- 
nels K a , indexed by a scale o\ More specifically, a e R + , 
J(:R 2 xR + ^ R + ; (x, o) ^K a (x). For convenience, we 
will consider the isotropic Gaussian kernel of the form 



Kr(x) 



I27Z 



exp 



'2a 2 



(2) 



We will work with the non-linear features proposed 
by Hong, Soatto and Vese [20], which were designed to 
retain boundary information. They are given by one of 
the following two expressions, 



R a :L\Q) -> L\R 2 ), 

S{x) h+ R a (x\S)=S(x)(K a * (1 - S{x))), 

or the symmetrized version 

R a : L 1 (to) -> L 1 (R 2 ) , 

S (x) R a (x |S) = S (x) {Ka * (1 - S (jc))) 

+ (l-S(x))(K ff *S(x)). 



(3) 



(4) 



The shape representation for different scales using the 
first (simpler) features, Eq. (3), are shown in Figure 3. 



The shape representation for different scales using the 
second (symmetric) features, Eq. (4), are shown in 
Figure 4. Both the binary representation (S) and the 
nonlinear shape features (7?^) include the original 
boundary information. However, the nonlinear shape 
features also encode the local shape information (up to 
a scale O"), which is not explicitly available in the binary 
representation. 

These shape features have several useful properties 
(for more details on these properties, please see [20]): (i) 
they are very robust under the presence of noise that 
gets incorporated in the segmentation process; (ii) since 
their values depend on the local geometry, these features 




Figure 4 Cardiocyte shape features (symmetric). Examples of the 
shape features R a (x \ S) = S (x) (K a * (1 - S (*))) + 0-5 (*)) (K a * S 
(x)) for (a) c= 15 and (b) o = 25. 
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propagate the shape information inside and outside the 
boundary; (iii) because the value of the features at a 
point is a local statistic of the shape in a neighborhood 
of that point, these features capture the context of the 
particular shape; and (iv) these shape features are very 
straightforward to compute. 

Cardiocyte Contractility Assessment via Shape Matching 

In this paper, the contractility analysis is done at the cel- 
lular level, thereby only measuring the overall contrac- 
tions in the cell. This is consistent with the contraction 
measurements that are being used in our laboratory. We 
are working on a similar energy conservation and infor- 
mation content approach for assessing contractility in 
neonatal cardiocytes, where we will measure the fine 
granular changes in the image. Unlike adult cardiocytes, 
which are highly organized and quite similar in morphol- 
ogy, the neonatal cardiocyte is in the process of develop- 
ing its contractile machinery. The neonatal cardiocyte is 
generally unable to retract its cell boundary during con- 
traction, and noticeable changes occur only within the 
cell perimeter. For these reasons, it is difficult to perform 
contractile measurements on this cell type in a manner 
similar to that of the adult cardiocyte, in which changes 
in cell boundary are quantified during contraction. 
(Please see a recent article by Bazan, Torres-Barba, Pao- 
lini and Blomgren [27] for a previously developed com- 
putational framework for the quantitative assessment of 
contractile responses of isolated neonatal cardiac 
myocytes.) 

We are given two shapes of the same topology, i.e., 
the shapes of a relaxed and contracted cardiocyte, 
respectively (Figure 5), defined by the two functions S lf 




Figure 5 Cardiocyte contraction/relaxation process. Example of 
two shapes of the same topology that depict the contraction/ 
relaxation process that we are trying to measure. 

v / 



S 2 : H — > {0, 1}. Our intent is to transform one into the 
other, and vice-versa, in a process that resembles that of 
the contraction/relaxation of the cardiocyte. As pro- 
posed in [20], we will do this by warping, that is a 
domain deformation h : O — > R 2 such that h (O) = O 
and 

Si (*) = S 2 (h(x)), Vx e Q. (5) 

We are interested in the distance between the two 
shapes, i.e., our proposed measure of contraction. The 
distance between the shapes will be defined as the 
energy of the aforementioned warping. Since there are 
infinitely many warpings that satisfy (5), in order to 
make this distance unique, it is defined as the one that 
minimizes the energy in a suitably chosen class [28]. For 
instance, as proposed in [20], d{S\, S2)=min/ 1 ||/z|| sub- 
ject to (5), where h is a diffeomorphism and ||-|| is 
some chosen norm of integral form. This minimization 
process can be recast as a variational problem of the 
form 

d{S lf S 2 ) = inf[E data (Si,S 2 |/i) +aE xeg (h)], (q 

heH v ' 

where H is a suitable function space. The distance 
between the two shapes, as defined above, is a function 
of two energy components: £data (Si, S 2 \h), which repre- 
sents the data fidelity; and £ reg (h), which is a regulari- 
zation term. Both terms are explained in detail below. 

In order to guarantee a meaningful vis-a-vis corre- 
spondence between the two shapes up to a certain scale 
[22] {i.e., a point x in the interior of the set that defines 
Si is mapped to a point h (x) in the interior of the set 
that defines S 2 ; and similarly, a point x in the exterior of 
the set that defines S x is mapped to a point h (x) in the 
exterior of the set that defines S 2 ), we adopt the mea- 
sure of data fidelity between the two shapes Si and S 2 
proposed in [20]: 

Edata(fr|Si,S 2 )= [ \Mx\Sl)-Ra{h{x)\S2)\ 2 dx t (7) 

where R G is the shape features (3) or (4). 

In our application, the purpose of the regularization 
term £ reg (h) is two-fold. First, it is there to render the 
problem well posed and is designed to penalize varia- 
tions of the diffeomorphism function h in favor of 
smoothness. Second, it makes the deformations compa- 
tible with the deforming material. Several authors have 
provided insight into the constitutive laws that describe 
the mechanical responses of resting and contracting car- 
diac muscle, along with their regional and temporal var- 
iations [23,29-31]. The problem is, however, extremely 
complex and has required of elaborate combinations of 
multiaxial tissue testing [32,33], microstructural 
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morphological modeling [34,35], statistical parameter 
estimation, and validation with measurements [36,37]. 
Very sophisticated numerical methods are also essential 
for accurate quantitative analysis in all phases of the 
investigations [23]. 

The intact cardiac muscle undergoes finite deforma- 
tions during the normal cardiac cycle. Thus, the classical 
linear theory of elasticity is inappropriate for resting 
myocardial mechanics [38,39,25]. The myocardium is 
frequently modeled as a finite hyper-elastic material, 
where the second Piola-Kirchhoff stress tensor compo- 
nents P t p are related to the components of the Lagran- 
gian Green's strain tensor E t p through the pseudo-strain 
energy W, as 



/3W dW\ , 



(8) 



where Cg is the right Cauchy-Green deformation tensor 
and p is a hydrostatic pressure Lagrange multiplier [23] 
(which we assume to be zero in this analysis). Several 
functional forms have been proposed for W [24,40,30, 
41,42,26,33]. We will adopt the transversely isotropic func- 
tional form proposed by Guccione, McCulloch and Wald- 
man [24] that considers the fibrous structure of the 
myocardium. The strain energy potential W is an expo- 
nential function of the strain components Ey referred to 
the fiber coordinates 



The optimal correspondence given by /z* is obtained by 
h* = argmin(£ da ta + E reg ). (12) 

The energy minimization is performed in a variational 
framework using a gradient descent method. The Euler- 
Lagrange equation corresponding to the energy E = 
£ data + £ reg yields the gradient direction for h: 



dh 

~dt 



dE 

dh 



dh 



BE, 



reg 



dh ' 



(13) 



where a is a small parameter (Lagrange multiplier). 
Using the features from Eq. (3), 
R a (x\S) = S (x) (K a • (1 - S (*))), we have 

= vs 2 °h ■ - MHx)\s 2 )) 

• {K 0 * {S 2 °h - 1)) (14) 
+ -Kir * {{Ra(x\Si) - Ra{h{x)\S 2 )) • S 2 °h)}, 

and, in components form, we have 



9£ re g 
~dhT 



-bfj^(2ui + 3u\ + 3u\ +v\+ v\ui) 



-b t -^(2u 2 v 2 + u 2 v\ + u\) 

— bf s j^(u\ + u 2 V\ + U\u\ + u 2 viv 2 ) 0-5) 

- bf s ^(u 2 + 2uiu 2 +vx+ viui + u\u 2 

+ V\V 2 + V\V 2 Ui)} r 



W= §(e Q - 1), 

Q = b f E 2 11+ b t E 2 22 + b fs {E 2 12+ E 2 21 ), 



(9) 



where 



1 /Shi dhj dh k dh k 

2 \dXj dXi dXi dXj 



(10) 



The aforementioned fiber coordinate system has the 
coordinate directions of the muscle fiber axis, the axis 
of the myofiber sheets, and the axis normal to the 
sheets, and is derived by rotating the cardiac coordinate 
system through the two angles that define the local 
myofiber-sheet orientation [24]. In Eq. (9), En is the 
fiber strain, E 22 is the cross-fiber in-plane strain, and E 12 
is the shear strain in the fiber-cross fiber coordinate 
plane. Omens, MacKenna and McCulloch [26] have 
found that the material constants C = 1.1 kPa, bf = 9.2 
kPa, b t = 2.0 kPa, and bf s = 3.7 kPa are appropriate to 
model the strains measured in the rat midwall. Then, 
the regularization term, £ reg (/z), can be written as 



1 Jq 



l)dx. 



(id 



— b t -^[2v 2 + 3v 2 +v\+u 2 2 + u 2 v 2 ) 

— bf s ^(u 2 + u 2 v 2 + v\ + 2v\V 2 + U\U 2 

+ U\U 2 V 2 +V\V 2 2 ) 

-bf s ^{viu 2 +v 2 x + V X U\U 2 + 1^V 2 )J , 

where 

Q = bf\{2u x +u 2 + v\) 2 + bt\[2v 2 +v\+ u\) 2 
+ bf s j(u 2 + v\ + U\U 2 + V\V 2 ) 2 , 

With the following short-hand notation 
(hlk = _ =Wl/ (hlk = _ =M2/ 

{h2)xi = — = Vll {h2)x2 = — = V2 . 



(16) 



(17) 



(18) 



Average Shape of the Relaxed State 

In order for us to use the energy of the warping as a 
measure of the cardiocyte's contractions, we need to 
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determine a baseline that will represent the state when 
the cell is not contracting (relaxed). We will define this 
baseline as the average shape of the cardiocyte's relaxed 
state. In other words, given an ensemble of shapes {Si, 
S 2 , S n }, we are interested in finding the average of the 
cardiocyte's shapes representing the uncontracted phase. 

A video from our lab depicting a contracting cardio- 
cyte will typically comprise about one thousand frames. 
On average, about ninety percent of these frames will 
show the cardiocyte in its relaxed state. 

Furthermore, the shapes of the cardiocyte in these 
frames are practically identical, then finding the average 
shape only guarantees a more unbiased measurement 
while providing for some regularization. Thus, in the 
interest of minimizing the overall computing time, we 
implemented a very simple averaging of contours in lieu 
of a shape averaging. The algorithm is as follows: 1) 
Identify the frames from each relaxed phase; 2) Obtain 
the contours of the shapes; 3) Resample the contours so 
that they will have the same number of points; 4) Find 
the average centroid point; 5) Interpolate the points in 
each contour using splines in polar coordinates; and 6) 
Interpolate the splines among the contours. Figure 6 
shows the effectiveness of this very simple contour aver- 
aging approach where the contours of a relaxed shape 
and a contracted shape were averaged. We show this 
averaging since the average of two uncontracted con- 
tours is practically indistinguishable from the two 
uncontracted contours. 

It is worth noting here that the aforementioned shape 
averaging could have been accomplished within the 
same shape deformation approach due to Hong, Soatto 
and Vese [20], where they define their shape average M 
as the shape that is closest to the ensemble, by simulta- 
neously minimizing their energy functional equivalent to 
Eq. (6). In other words, they look for M that minimizes 

n 

J>(M,S,). (19) 

1=1 




Figure 6 Average relaxed shape. The contours of a relaxed shape 
and a contracted shape were averaged in order to show the 
effectiveness of the contour averaging method. In practice, we only 
average the contours of uncontracted shapes. 

V / 



They perform this optimization via alternating mini- 
mization and gradient descent. Implementing the opti- 
mization process proposed in [20] was deemed to be 
too expensive a proposition for our purposes. Thus, we 
opted for performing the simpler average of contours 
which is two orders of magnitudes less expensive to run 
on a personal computer equipped with MATLAB™. 

Results and Discussion 

Experimental Results 

A typical deformation of the shape of a contracting myo- 
cyte into the shape of the average relaxed myocyte is 
shown in Figure 7. In this particular example, the warp- 
ing process follows the path of minimum energy (as 
defined in the Methods section) and the warped shape 
matches the average shape with a correlation coefficient 
of over 99.99%. The perfect matching was obtained in 
246 iterations after which the warping energy reached a 
stable minimum. This process is repeated for every frame 
in the video depicting the contracting myocyte and the 
total energy employed in each deformation is calculated. 
These energies represent our measure of contractility. 




o 1 1 1 1 1 1 

0 50 100 150 200 250 

Iterations 



Figure 7 Warping process and warping energy (a) Warping of 

the shape of a contracting myocyte into the shape of the average 

relaxed myocyte. After 246 iterations the template shape matched 

the target shape with a similarity of over 99.99% as measured by 

the correlation coefficient between the warped shape and the 

average shape, (b) Energy of the warping of the shape of a 

contracting myocyte into the shape of the average relaxed 

myocyte. The energy minimization process reaches a stable 

minimum after 246 iterations, 
k J 
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Figure 8 shows the normalized contraction measures 
along with the energy profile resulted from the deforma- 
tion process. 

We tested the proposed approach by assessing the con- 
tractile responses in isolated adult rat cardiocytes har- 
vested and imaged as described in [18]. We contrasted 
these measurements against both the classic raster-line 
approach [15,11,12] and the contractility pipeline 
described by Bazan, Torres-Barba, Paolini and Blomgren 
[18]. Our results show good qualitative and quantitative 
agreements between the proposed method and both the 
raster-line method and the contractility pipeline as far as 
frequency, pacing, and overall behavior of the contrac- 
tions are concerned. Figure 9 reproduces the average 
normalized contractions as assessed by the three meth- 
ods in this study. We observe great similarities among 
the three methods, specially during the contraction or 
activation process marked by the electrical stimulus. 

There exist small discrepancies during the relaxation 
phase where the raster-line method seems to show a 



slightly different recover process. This phenomenon was 
already reported in [18]. The raster-line method-being a 
one-dimensional technique-is unable to capture the full 
extent of the contraction process occurring outside its 
domain of influence. The proposed method, as well as 
the computational pipeline, capture the contraction of 
the cell as a two dimensional event over the entire 
boundary of the cell. Imaging the contracting cells with 
a high speed camera might eventually elucidate this 
small disagreement. 

Conclusions 

We explored the viability of a new approach to cardiocyte 
contractility assessment based on biomechanical proper- 
ties of the cardiac cells, energy conservation principles, 
and information content measures. We defined our mea- 
sure of cell's contraction as being the distance between 
the shapes of the contracting cell, assessed by the total 
energy of the domain deformation (warping) of one cell 
shape into another. To guarantee a meaningful vis-a-vis 
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Figure 9 Average Contraction. Average contraction (10 contractions 
in total) as assessed by the three methods in this study: proposed 
method (red line), computational pipeline (green line), and raster-line 
(blue line). 



correspondence between the two shapes, we employed 
both a data fidelity term and a regularization term. The 
data fidelity term is based on nonlinear features of the 
shapes while the regularization term enforces the com- 
patibility between the shape deformations and that of a 
hyper-elastic material Our results show good qualitative 
and quantitative agreements between the proposed 
method and both the raster-line method and the contrac- 
tility pipeline as far as frequency, pacing, and overall 
behavior of the contractions are concerned. 

We hypothesize that this methodology, once appropri- 
ately developed and customized, can provide a frame- 
work for computational cardiac cell bio-mechanics that 
can be used to integrate both theory and experiment. For 
example, besides giving a good assessment of contractile 
response of the cardiocyte, since the excitation process of 
the cell is a closed system, this methodology can be used 
in an attempt to infer statistically significant model para- 
meters for the constitutive equations of the cardiocytes. 
This conjecture is still very preliminary. Nonetheless, the 
way we envision this analysis resembles that of finding 
the "spring constant" by measuring the deformation in 
the material under a constant load. In our case, we know 
the deformation of the cell undergoing contraction as 
measured by the changes in the shape. We also know 
that the energy that gets incorporated into the system 
through the electrical stimulus is the same for every con- 
traction. For the sake of this argument, let us assume 
that it takes 10 snap-shots for the cell to go from its 
relaxed state to the fully contracted state (that is 10 
shapes S 0 , Si, S 8 , <S 9 ). Then, the total energy employed 
by the cell to go from S 0 to S 9 has to equal the sum of 
the energies for moving the cell from S 0 to Si + Si to S 2 



+ ...+ S 7 to S 8 + S 8 to S 9 . If we do this for every contrac- 
tion in the experiment, we will have an over-determined 
system that we can solve (one of the solutions) by finding 
the least-square error. Note that, since we are assuming 
that the constitutive parameters are the same for the tis- 
sue sample, we can simply average these parameters 
across many cells in the same experiment. 

The aforementioned possible integration between the- 
ory and experiment can also be extended further to 
include functional coupling between the many physiolo- 
gical processes that interact with mechanics such as cell 
growth and signaling, metabolism, transport and 
electrophysiology. 

Acknowledgements 

Carlos Bazan would like to thank Dr. Byung-Woo Hong and Dr. Andrew 
McCulloch for their advice. The authors thank Xian Zhang for her help in the 
preparation of biological procedures. This work has been supported in part 
by NIH Roadmap Initiative award R90 DK07015 and NIH NIDDK, the 
California Metabolic Research Foundation, and the Computational Science 
Research Center at San Diego Sate University. 

Author details 

Computational Science Research Center, San Diego State University, 5500 
Campanile Drive, San Diego, CA 92182-1245, USA. department of 
Mathematics & Statistics, San Diego State University, 5500 Campanile Drive, 
San Diego, CA 92182-7720, USA. department of Biology, San Diego State 
University, 5500 Campanile Drive, San Diego, CA 92182-4614, USA. 

Authors' contributions 

CB conceived the idea and designed the methodology and the experiments; 
DTB conducted the biological data acquisition; PB supplied key steps for the 
calculus of variations-then executed by CB and verified by TH; CB and TH 
conducted the numerical computations; PP provided biological insight and 
veracity to the project; CB wrote the manuscript; All authors read and 
approved the final manuscript. 

Received: 1 1 June 201 1 Accepted: 22 August 201 1 
Published: 22 August 2011 

References 

1. Banyasz T, Lozinskiy H, Payne C, S E, Norton B, Chen B, Y CI, Izu L, Balke C: 
Transformation of adult rat cardiac myocytes in primary culture. 
Experimental Physiology 2007, 93(3):370-382. 

2. Del bridge L, Roos K: Optical methods to evaluate the contractile function of 
unloaded isolated cardiac myocytes. Molecular Cell Cardiology 1 997, 29:1 1 -1 5. 

3. Maltsev V, HN S, Tanimura M, Lesch M, Goldstein S, Undrovinas A: 
Relationship between action potential, contraction-relaxation pattern, 
and intracellular Ca2+ transient in cardiomyocytes of dogs with chronic 
heart failure. Cell Mol Life Sciences 1998, 54:597-605. 

4. Mitcheson JS, Hancox J, Levi A: Cultured adult cardiac myocytes: future 
applications, culture methods, morphological and electrophysiological 
properties. Cardiovascular Research 1 998, 39:280-300. 

5. Ren J, Wold L: Measurement of cardiac mechanical function in isolated 
ventricular myocytes from rats and mice by computerized video-based 
imaging. Biology Procedures Online 2001, 3. 

6. Shevchuk A, Gorelik J, Harding S, Lab M, Klenerman D, Korchev Y: 
Simultaneous measurement of Ca 2+ and cellular dynamics: Combined 
scanning ion conductance and optical microscopy to study contracting 
cardiac myocytes. Biophysical Journal 2001, 81:1759-1764. 

7. Lecarpentier Y, Martin J, Claes V, Chambaret J, Migus A, Antonetti A, Hatt P: 
Real-time kinetics of sarcomere relaxation by laser diffraction. Circulation 
Research 1985, 56:331-339. 

8. Philips C, Duthinh Y, Houser S: A simple technique to measure the rate 
and magnitude of shortening of single isolated cardiac myocytes. IEEE 
Transactions on Biomedical Engineering 1986, 33:929-934. 



Bazan et al. BMC Biophysics 201 1, 4:1 7 Page 9 of 9 

http://www.biomedcentral.com/2046-1 682/4/1 7 



10. 



12. 



14. 



15. 



16. 



17. 



19. 



20. 



21. 



22. 



23. 



24. 



25. 



26. 



27. 



28. 
29. 



30. 



32. 



33. 



34. 



Mukherjee R, Crawford F, Hewett K, Spinale F: Cell and sarcomere 
contractile performance from the same cardiocyte using video 
microscopy. Journal of Applied Physiology 1993, 74:2023-2033. 
Rieser G, Sabbadini R, Paolini P, Fry GM: Inesi: Sarcomere motion in 
isolated cardiac cells. AJP - Cell Physiology 1979, 236:C70-C77. 
Steadman B, Moore K, Spitzer K, Bridge J: A video system for measuring 
motion in contracting heart cells. IEEE Transactions on Biomedical 
Engineering 1988, 35(4):264-272. 

Tameyasu T, Toyoki T, Sugi H: Nonsteady motion in unloaded 
contractions of single frog cardiac cells. Biophysical Journal 1985, 
48:461-465. 

Krueger J, Denton A, Siciliano G: Nature of motions between sarcomeres 
in asynchronously contracting cardiac muscle cells. Biophysical Journal 
1992, 61:145-160. 

Paolini P, Roos K, Baskin R: Light diffraction studies of sarcomere 
dynamics in single skeletal muscle fibers. Biophysical Journal 1977, 

20(2):221-232. 

Harris P, Stewart D, Cullinan M, Delbridge L, Dally L, Grinwald P: Rapid 
measurement of isolated cardiac muscle cell length using a line-scan 
camera. IEEE Transactions on Biomedical Engineering 1987, 34:463-467. 
Murphy D: Fundamentals of light microscopy and electronic imaging New 
York, New York: Wiley-Liss; 2001. 

Reuben H, Godinez M, Paolini P, Bejar E: Analysis of contractile dynamics 
of adult cardiac myocytes using a computer controlled edge detector: 
Effects of 5-hydroxytryptamine. Cardiovascular Pathobiology 1998, 
2:149-158. 

Bazan C, Torres D, Paolini P, Blomgren P: Image processing techniques for 
assessing contractility in isolated adult cardiac myocytes. International 
Journal of Biomedical Imaging 2009, 2009(52954):1 1. 
Mukherjee R, Spinale F, Crawford F: Measurement of dynamic cellular and 
sarcomere contractile properties from the same cardiocyte. Proceedings 
of Annual International Conference of the IEEE Engineering in Medicine and 
Biology Society 1992, 2:392-393. 

Hong BW, Soatto S, Vese L: Enforcing local context into shape statistics. 

Advances in Computational Mathematics 2009, 31:185-213. 
Hong BW, Prados E, Vese L, Soatto S: Shape representation based on 
integral kernels: Application to image matching and segmentation. IEEE 
Conference on Computer Vision and Pattern Recognition, IEEE 2006, 833-840. 
Manay S, Cremers D, Hong BW, Yezzi A, Soatto S: Integral invariants for 
shape matching. IEEE Transactions on Pattern Analysis and Machine 
Intelligence 2006, 28(1 0):1 602-1 61 8. 

Costa K, Holmes J, McCulloch A: Modeling cardiac mechanical properties 

in three dimensions. Philosophical Transactions of the Royal Society London 
A 2001, 359:1233-1250. 

Guccione J, McCulloch A, Waldman L: Passive material properties of intact 
ventricular myocardium determined from a cylindrical model. Journal of 
Biomechanical Engineering 1991, 1 18(2):262-263. 

McCulloch A: In Cardiac Biomechanics, Biomedical Engineering Fundamentals. 
Volume chap 54. San Diego, California: CRC Press; 2006. 
Omens J, MacKenna D, McCulloch A: Measurement of strain and analysis 
of stress in resting rat left ventricular myocardium. Journal of 
Biomechanics 1993, 26(6):665-676. 

Bazan C, Torres D, Paolini P, Blomgren P: Image processing techniques for 
the assessment of contractile responses in neonatal cardiac myocytes. 

International Journal of Biomedical Imaging 201 1, 201 1:9, (Article ID 729732). 
Grenander U: General Pattern Theory Oxford: Oxford University Press; 1993. 
Demer L, Yin F: Passive biaxial mechanical properties of isolated canine 
myocardium. Journal of Physiology 1983, 339:615-630. 
Humphrey J, Strumpf R, Yin F: Determination of a constitutive relation for 
passive myocardium: II. Parameter estimation. Journal of Biomechanical 
Engineering 1 990, 1 1 2(3):340-346. 

LeGrice I, Hunter P, Young A, Smail B: The architecture of the heart: A 

database model. Philosophical Transactions of the Royal Society London A 
2001, 359:1217-1232. 

Sacks M, Sun W: Multiaxial mechanical behavior of biological materials. 

Annual Review of Biomedical Engineering 2003, 5:251284. 

Yin F, Strumpf R, Chew P, Zeger S: Quantification of the mechanical 

properties of non contracting canine myocardium under simultaneous 

biaxial loading. Journal of Biomechanics 1987, 20(6):577589. 

Jacquemet V, Henriquez C: Loading effect of broblast-myocyte coupling 

on resting potential, impulse propagation, and repolarization: Insights 



35. 



36. 



37. 



40. 



42. 



from a microstructure model. American Journal of Physiology - Heart and 
Circulatory Physiology 2008, 294(5):H2040-H2052. 
Roberts S: Non-uniform interstitial loading in cardiac microstructure 
during impulse propagation. PhD thesis Duke University, Durham, North 
Carolina; 2009. 

Bursac N, Papadaki M, Cohen R, Schoen F, Eisenberg S, Carrier R, Vunjak- 
Novakovic G, Freed L: Cardiac muscle tissue engineering: Toward an in 
vitro model for electrophysiological studies. American Journal of 
Physiology - Heart and Circulatory Physiology 1 999, 277(2):H433-H444. 
Vliegen H, Van Der Laarse A, Cornelisse C, Eulderink F: Myocardial changes 
in pressure overload-induced left ventricular hypertrophy. A study on 
tissue composition, polyploidization and multinucleation. European Heart 
Journal 1991, 1 2(4):488-494. 

Cowin S, Doty S: Tissue Mechanics New York, New York: Springer; 2007. 
Fung Y: Biomechanics: Mechanical Properties of Living Tissues. 2 edition. 
Germany: Springer-Verlag, Inc; 1993. 

Gupta K, Ratcliffe M, Fallert M, Edmunds L Jr, Bogen D: Changes in passive 
mechanical stiffness of myocardial tissue with aneurysm formation. 

Circulation 1994, 89:2315-2326. 

Humphrey J, Yin F: A new constitutive formulation for characterizing the 
mechanical behavior of soft tissues. Biophysical Journal 1987, 
52(4):563-570. 

Novak V, Yin F, Humphrey J: Regional mechanical properties of passive 

myocardium. Journal of Biomechanics 1994, 27(4):403-412. 



doi:1 0.1 1 86/2046-1 682-4-1 7 

Cite this article as: Bazan et al:. Introduction of non-linear elasticity 
models for characterization of shape and deformation statistics: 
application to contractility assessment of isolated adult cardiocytes. 

BMC Biophysics 2011 4:17. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at 
www.biomedcentral.com/submit 



o 



BioMed Central 



